Search for: All records

Abstract The first infall of the LMC into the Milky Way (MW) represents a large and recent disruption to the MW circumgalactic medium (CGM). In this work, we use idealized, hydrodynamical simulations of an MW-like CGM embedded in a dark matter halo with an infalling LMC-like satellite initialized with its own CGM to understand how the encounter is shaping the global physical and kinematic properties of the MW CGM. First, we find that the LMC drives order-unity enhancements in MW CGM density, temperature, and pressure due to a $\mathcal{M}\approx 2$shock from the supersonic CGM–CGM collision. The resulting shock front extends from the LMC to beyond ∼R_200,MW, amplifying column densities, X-ray brightness, thermal Sunyaev–Zeldovich distortion, and potentially synchrotron emission from cosmic rays over large angular scales across the southern hemisphere. Second, the MW’s reflex motion relative to its outer halo induces a dipole in CGM radial velocities, withv_R ± 30–50 km s⁻¹atR > 50 kpc in the northern and southern hemispheres, respectively, consistent with measurements in the stellar halo. Finally, ram pressure strips most of the LMC’s CGM, leaving ∼10⁸⁻⁹M_⊙warm ionized gas along the past orbit of the LMC, moving at high radial and/or tangential velocities ∼50–100 kpc from the MW. Massive satellites like the LMC leave their mark on the CGM structure of their host galaxies, and signatures of such interactions may be observable in key all-sky tracers of the MW CGM and those of other massive galaxies. 

Abstract The large-scale morphology of Milky Way (MW)–mass dark matter halos is shaped by two key processes: filamentary accretion from the cosmic web and interactions with massive satellites. Disentangling their contributions is essential for understanding galaxy evolution and constructing accurate mass models of the MW. We analyze the time-dependent structure of MW-mass halos from zoomed cosmological-hydrodynamical simulations by decomposing their mass distribution into spherical harmonic expansions. We find that the dipole and quadrupole moments dominate the gravitational power spectrum, encoding key information about the halo’s shape and its interaction with the cosmic environment. While the dipole reflects transient perturbations from infalling satellites and damps on dynamical timescales, the quadrupole—linked to the halo’s triaxiality—is a persistent feature. We show that the quadrupole’s orientation aligns with the largest filaments, imprinting a long-lived memory on the halo’s morphology even in its inner regions (∼30 kpc). At the virial radius, the quadrupole distortion can reach 1–2 times the spherical density, highlighting the importance of environment in shaping MW-mass halos. Using multichannel singular spectrum analysis, we successfully disentangle the effects of satellite mergers and filamentary accretion on quadrupole. We find that, compared to isolated MW–LMC simulations that typically use a spherical halo, the LMC-mass satellite induces a quadrupolar response that is an order of magnitude larger in our cosmological halo. This highlights the need for models that incorporate the MW’s asymmetry and time evolution, with direct consequences for observable structures such as disk warps, the LMC-induced wake, and stellar tracers—particularly in the era of precision astrometry. 

Abstract The infall of the Large Magellanic Cloud (LMC) into the Milky Way’s halo impacts the distribution of stars and dark matter (DM) in our Galaxy. Mapping the observational consequences of this encounter can inform us about the properties of both galaxies, details of their interaction, and possibly distinguish between different DM models.N-body simulations predict a localized overdensity trailing the LMC’s orbit both in baryonic and DM, known as the wake. We collected wide-field, deep near-infrared, and optical photometry using VIRCAM and DECam across four fields along the expected wake, covering the sky region expected to span most of its predicted density contrast. We identify over 400 stars comprising two different tracers, near main-sequence turnoff stars and red giants, which map the halo between 60 and 100 kpc, deriving stellar halo densities as a function of sky position and Galactocentric radius. We detect (1) a break in the halo radial density profile at 70 kpc not seen in northern halo studies and (2) a clear halo overdensity starting also at 70 kpc, with density contrast increasing steadily toward the expected current location of the wake. If this overdensity is the LMC wake, its peak density contrast is as pronounced as the most massive LMC model considered. Contamination from unidentified substructures may bias our wake detections, so wider-area surveys with similar depth are needed for confirmation. 

We construct time-evolving gravitational potential models for a Milky Way–mass galaxy from the FIRE-2 suite of cosmological-baryonic simulations using basis function expansions. These models capture the angular variation with spherical harmonics for the halo and azimuthal harmonics for the disk, and the radial or meridional plane variation with splines. We fit low-order expansions (four angular/harmonic terms) to the galaxy’s potential for each snapshot, spaced roughly 25 Myr apart, over the last 4 Gyr of its evolution, then extract the forces at discrete times and interpolate them between adjacent snapshots for forward orbit integration. Our method reconstructs the forces felt by simulation particles with high fidelity, with 95% of both stars and dark matter, outside of self-gravitating subhalos, exhibiting errors ≤4% in both the disk and the halo. Imposing symmetry on the model systematically increases these errors, particularly for disk particles, which show greater sensitivity to imposed symmetries. The majority of orbits recovered using the models exhibit positional errors ≤10% for 2–3 orbital periods, with higher errors for orbits that spend more time near the galactic center. Approximate integrals of motion are retrieved with high accuracy even with a larger potential sampling interval of 200 Myr. After 4 Gyr of integration, 43% and 70% of orbits have total energy and angular momentum errors within 10%, respectively. Consequently, there is higher reliability in orbital shape parameters such as pericenters and apocenters, with errors ∼10% even after multiple orbital periods. These techniques have diverse applications, including studying satellite disruption in cosmological contexts. 

Abstract Stellar streams in the Milky Way are promising detectors of low-mass dark matter (DM) subhalos predicted by ΛCDM. Passing subhalos induce perturbations in streams that indicate the presence of the subhalos. Understanding how known DM-dominated satellites impact streams is a crucial step toward using stream perturbations to constrain the properties of dark perturbers. Here, we cross-match a Gaia Early Data Release 3 and SEGUE member catalog of the Cetus-Palca stream (CPS) with H3 for additional radial velocity measurements and fit the orbit of the CPS using this six-dimensional (6D) data. We demonstrate for the first time that the ultra-faint dwarf Segue 2 had a recent (77 ± 5 Myr ago) close flyby (within the stream's 2σwidth) with the CPS. This interaction enables constraints on Segue 2’s mass and density profile at larger radii ( $O\left(1\right)$kpc) than are probed by its stars ( $O\left(10\right)$pc). While Segue 2 is not expected to strongly affect the portion of the stream covered by our 6D data, we predict that if Segue 2’s mass within ∼ 6 kpc is 5 × 10⁹M_⊙, the CPS's velocity dispersion will be ∼ 40 km s⁻¹larger atϕ₁ > 20° than atϕ₁ < 0°. If no such heating is detected, Segue 2’s mass cannot exceed 10⁹M_⊙within ∼ 6 kpc. The proper motion distribution of the CPS near the impact site is mildly sensitive to the shape of Segue 2’s density profile. This study presents a critical test for frameworks designed to constrain properties of dark subhalos from stream perturbations. 

Abstract We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulations of the MW-LMC system, we find that the DM particles in the Solar neighborhood originating from the LMC analogue dominate the high speed tail of the local DM speed distribution. Furthermore, the native DM particles of the MW in the Solar region are boosted to higher speeds as a result of a response to the LMC's motion.We simulate the signals expected in near future xenon, germanium, and silicon direct detection experiments, considering DM interactions with target nuclei or electrons. We find that the presence of the LMC causes a considerable shift in the expected direct detection exclusion limits towards smaller cross sections and DM masses, with the effect being more prominent for low mass DM. Hence, our study shows, for the first time, that the LMC's influence on the local DM distribution is significant even in fully cosmological MW analogues. 

ABSTRACT Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e. proto-galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: (i) can have a stellar mass range between 1 × 108 < M⋆ < 2 × 1010 [M⊙], a virial mass range between 3 × 1010 < M⋆ < 6 × 1011 [M⊙], and be as young as 8 ≲ Age ≲ 12.8 [Gyr] (1 ≲ z ≲ 6); (ii) are pre-dominantly centrally concentrated, with $$\sim 50~{{\ \rm per\ cent}}$$ of the stars contained within 5–10 kpc; (iii) on average show weak but systematic net rotation in the plane of the host’s disc at z = 0 (i.e. 0.25 ≲ 〈κ/κdisc〉 ≲ 0.8); (iv) present [α/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host’s old disc; and (v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that $$\sim 60~{{\ \rm per\ cent}}$$ of the proto-Milky Way galaxies are comprised by 1 dominant system (1/5 ≲M⋆/M⋆, proto-MilkyWay≲ 4/5) and 4–5 lower mass systems (M⋆/M⋆, proto-MilkyWay≲ 1/10); the other $$\sim 40~{{\ \rm per\ cent}}$$ are comprised by 2 dominant systems and 3–4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems. 

Abstract The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW’s stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, versus fuzzy DM (FDM) withm_a= 10⁻²³eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake’s density by ∼10%, and holds the wake together over larger distances (∼50 kpc) than if self-gravity is ignored. The DM wake’s mass is comparable to the LMC’s infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is ∼20% dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars’ kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM versus FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70–100 kpc. 

Abstract The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the “Timing Argument,” which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history and including the recent pericentric passage of the Large Magellanic Cloud (LMC), and recent work has found that the MW disk is moving with a lower bound “travel velocity” of ∼32 km s⁻¹with respect to the outer stellar halo. Previous Timing Argument measurements have attempted to account for this nonequilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10%–12% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume. 

Abstract In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological-baryonic simulations. Using these potentials, we calculate actions for star particles in tidal streams around three galaxies with varying merger histories at each snapshot from 7 Gyr ago to the present day. We determine the action-space coherence preserved by each model using the Kullback–Leibler divergence to gauge the degree of clustering in actions and the relative stability of the clusters over time. We find that all models produce a clustered action space for simulations with no significant mergers. However, a massive (mass ratio prior to infall more similar than 1:8) interacting galaxy not present in the model will result in mischaracterized orbits for stars most affected by the interaction. The locations of the action space clusters (i.e., the orbits of the stream stars) are only preserved by the time-evolving model, while the time-independent models can lose significant amounts of information as soon as 0.5–1 Gyr ago, even if the system does not undergo a significant merger. Our results imply that reverse-integration of stream orbits in the MW using a fixed potential is likely to give incorrect results if integrated longer than 0.5 Gyr into the past.